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ABSTRACT 

Recent data support the idea that the filaments observed in Ha emission near the 
centres of some galaxy clusters were shaped by bulk flows within their intracluster 
media. We present numerical simulations of evaporated clump material interacting 
with impinging winds to investigate this possibility. In each simulation, a clump falls 
due to gravity while the drag of a wind retards the fall of evaporated material leading 
to elongation of the tail. However, we find that long filaments can only form if the 
outflowing wind velocity is sufficiently large, ~ 10*cms~^. Otherwise, the tail ma- 
terial sinks almost as quickly as the cloud. For reasonable values of parameters, the 
morphological structure of a tail is qualitatively similar to those observed in clusters. 
Under certain conditions, the kinematics of the tail resemble those reported in Hatch 
et al.(2006). A comparison of the observations with the numerical results indicates 
that the filaments are likely to be a few tens of Myrs old. We also present arguments 
which suggest that the momentum transfer, from an outflowing wind, in the forma- 
tion of these filaments is probably significant. As a result, tail formation could play a 
role in dissipating some of the energy injected by a central AGN close to the cluster 
centre where it is needed most. The trapping of energy by the cold gas may provide 
an additional feedback mechanism that helps to regulate the heating of the central 
regions of galaxy clusters and couple the AGN to the ICM. 

Key virords: 



1 INTRODUCTION 

Optical emission-line nebulae commonly surround massive 
alaxies in the c entre s of X-ray bright, cool cluster cores 
Crawford et al ] 119991 ). The origin of the Ha filaments 
has been attributed to a variety of processes including: 
condensation from an intracluster medium (ICM) evolv- 
ing as a cooling flow (iFabian et al. 1 1 1984 iHeckman et al.l 



1 19891 : iDonahue fc Voit 199 j): accretio n of clouds captured 
in galaxy merg ers (iBraine et al.l 1995); exp ulsion from the 
central galaxy (|Burbidge fc Burbidgel 19651 ) . 

NGC 1275 at the centre of the Perseus cluster con- 
tains the best studied example of such a nebula (|Hatch et al.l 
I2OO6I ). Its filaments are typically 50-100 pc thick and up to 30 
kpc lo ng, and the majority of them are radial. iHatch et al.l 
l|2006l) presented kinematic data that rules out dynamical 
models of purely infalling filaments. The observed kinematic 
properties provide strong evidence that the filaments are not 
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in gravitational free fall, because, if they were, their veloc- 
ities would rise sharp ly towards the centre of the nebula 
l|Heckman et al.l Il989l ) . The most conclusive evidence lies 
in the velocity structures of the northern and northwestern 
filaments. The lower half of the northern filament is red- 
shifted with respect to the galaxy, whilst the upper section 
is blueshifted. Thus, the upper part of the filament is flow- 
ing away from the galaxy while the lower part is flowing into 
the galaxy. 

The filaments may follow streamlines of the fiow of more 
tenuous material around them. In the Perseus cluster, some 
filaments appear to have been shaped by the wakes of the 
buoyant bubbles. Given that the data suggest that, at least, 
some parts of the filaments are outflowing, their origin may 
lie within the galaxy. NGC 1275 contains a large reservoir 
of cold molecular gas that could fuel them. 

The filaments have morphologies similar to those of 
laminar jets, and their relatively smooth structures may 
place constraints on turbulent motions in the ICM. Alter- 
natively, an ordered, amplified magnetic field trailing be- 
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hind a buoyant bubble interacting with filaments may pre- 
vent the destruction of the filaments by turbulent motions 
(|Ruszkowski et alJlioOTh . 

The optical nebula around NGC 1275 emits 4.1 x 
lO^^ergs"^ in Ha and N[II]. The nature of the power source 
remains unclear. Various excitation mechanisms for these 
lines have been proposed. For instance, although ionisation 
by the central AGN may be important for the inner re- 
gions, it cannot be the dominant source of power for the 
extended nebula because the H a intensity does not decreas e 
with distance from the nucleus l| Johnstone fc Fabianlll98^ . 
The nebula may be excited by stellar UV, but there is no 
spatial correlation between the filaments and the stellar clus- 
ters. Excitation by X-rays from the ICM seems unlikely, as 
the ICM may be as much as a hundred times le ss luminous 
in UV than in the X-rays (|Fabian et aLlbOOj ). Heat con- 
duction fr om the ICM to the co lder filaments has also been 
proposed (iDonahue et aklboOOl ). but it might also lead to 
the evaporation of the filaments on too short a timescale 
l|Nipoti fc Binnevll2004l ) . Shocks and turbulent m ixing lay- 
ers might play a role (|Crawford fc FabianI [l993 ) . Cosmic 
rays, preferentially diffusing along the magnetic field lines 
trailing behind rising bubbles, could possibly drive the ex- 
citation in those fi l amen ts that are located in bubble wakes 
l|Ruszkowski et al] |2007| ). The same magnetic fields lines 
could also prevent the filaments from evaporating due to 
thermal conduction. 

The same mechanisms that power the emission could 
also have a profound effect on the morphology of the emit- 
ting region. It is impossible to take all of these processes into 
account. Consequently, we will investigate the effect of grav- 
ity and an impinging wind that strips material from a cloud 
on the morphology and kinematics of the resulting filament. 

The main aim of the work reported here is the calcula- 
tion of the density and velocity structures of material evap- 
orated from clumps embedded in winds from the central 
galaxies of galaxy clusters and the comparison of model re- 
sults with observations. Section 2 contains some preliminary 
considerations, while the model assumptions and numerical 
approach are treated in section 3. The simulation results 
appear in section 4. In section 5 we investigate the devel- 
opment of tails in realistic environments by including the 
appropriate density and gravity for a galaxy cluster. Sec- 
tion 6 concerns momentum tranfer between clumps and the 
winds surrounding them. Section 7 discusses the possible 
trends in optical emission between clusters, and section 8 is 
a summary. 



2 PRELIMINARY CONSIDERATIONS 
2.1 Filamentary widths 

We consider a clump, or cloud, embedded in a lower-density, 
hotter ambient wind with a speed in the frame of the 
central galaxy. Material evaporates from the clump and in- 
teracts with the wind to form a filament. 

It may be possible to estimate the age of an observed 
filament from its width. Previous numerical studies have 
shown that for a broad range of parameters the width of 
a filament is comparable to the diameter of the clu mp from 
which material evaporates (e.g. iDvson et al]|2006l ). The ra- 



dius of the filament rf would have a lower bound of approx- 
imately: rf ~ 0.5cs,cAt, where At is the age of the filament, 
and Cs,c is the sound speed in the clump at the head of the fil- 
ament. The length of the filament U would have a maximum 
given roughly by l{ ~ 2rf«w/cs,c. The largest observed h/rt 
ratios would require that Vw/cs.c ~ 100s which is plausible 
if the clump material is molecular (and, thus, cold) before 
being heated to the filament temperature. 

Alternatively, the width of the filament may be gov- 
erned by the diffusion of material into the ICM. Conse- 
quently, it is possible to set an upper bound to the diffu- 
sion coefficient. The root mean square (RMS) displacement 
of a particle after time t is yj (r^) = V6Dt, where D is the 
diffusion coefficient. If the RMS displacement due to diffu- 
sion equals the measured radius of a filament, the filament's 
length is approximately 



bV 



(1) 



Substituting order of magnitude estimates: k = 30 kpc, 
rf = 100 pc, Uw = lO^kms"^ we find that D ~ lO^'^cm^ s"^ 
Interestin gly, this is comparable to the diffusion coefficients 
found by Roediger et al.l l|2007l ) at 30 kpc in their simula- 



tions of metal diffusion in the Perseus cluster. If diffusion 
were contributing significantly to filament broadening, this 
value of the diffusion coefficient would require the filaments 
to be roughly At ~ 5 x 10^ yrs old. 

This diffusion coefficient greatly exceeds the Spitzer 
thermal and viscous diffusivities. Turbulent diffusion would 
have to be important for such a large diffusion coefficient to 
obtain. The turbulent diffusion coefficient is given approx- 
imately by D ~ vl/3 where v is the characteristic veloc- 
ity and I the characteristic length scale of the largest ed- 
dies. If the turbulence is driven by the shear of the fluid 
flow past the clump and along the filament, we might as- 
sume that the thickness of the boun dary layer is roughly 3 
% of the length of t he filament (jHartguist fc DvsonlllQSa : 
ICanto fc Ragalll99ll ) which exceeds the observed filament 
widths. Therefore, a diffusion picture is not appropriate for 
the dynamics of the filament, but clearly viscous coupling 
across the entire cross section of a filament is strong, lead- 
ing to a fiow velocity that depends on the distance along the 
filament but not on the location relative to the nearest edge. 



2.2 Clump mass-loss 

The rate at which mass is ablated from a clump, m, might be 
estimated by the rate at which momentum associated with 
the wind is transferred to the clump, divided by a speed, Cs,a, 
anywhere between the sound speed of material in the clump 
and the difference between the average clump speed and the 
average speed of material in the filament. The appropriate 
value of Cs,a depends on how effective the viscous coupling 
(due to turbulence or any other mechanism) between the 
filament material and the wind material is. The total mass 
lost by a clump would then be given by 



rhAt = At 



(2) 



where A is the clump cross section and pw is the mass den- 
sity of the wind. If At is 30 Myr, the clump has a radius 
fc ~ 50 pc, pw ~ lO'^'^g cm~^, «w ~ lO'^km s~^ and 
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Cs,a ~ 0.3km s~^, the above equation gives m ^ 10* Mq. 
This is comparable to the mass of a structure with a num- 
ber density of 10 cm~'^, radius of 50 pc, and length of 30 
kpc. 



either case, an outflowing wind would ensure both infalling 
and outflowing portions of a filament. 



2.3 Clump motion 

A model of the clump's motion taking into account the loss 
of mass and the ram pressure of the impinging wind gives, 

d 2 
— {mx) = —mx — mg + KAp„{vvj~x) . (3) 

K is £L multiplicative factor which in the simplest case would 
be a constant of order unity, g is the gravitational field 
strength, x is the height of the clump above the central 
galaxy's mid-plane, and a dot indicates differentiation with 
respect to time. Expansion of the left-hand-side of equation 
((3| demonstrates that mass loss from the clump does not 
affect the motion of the clump. Therefore, the terminal ve- 
locity of the clump is. 

Using the values given above for pw and Vc and taking 
g ~ 10"*cm s"^, K = 1, and m ~ IO^Mq, we find that 
{mg/ KApwY''^ ~ lO'^kms^^, and so is significant for slower 
outfiow velocities. For larger masses, and gravities, the ter- 
minal velocity is greater. The terminal velocity of the clump 
is not likely to be constant, because the mass of the clump 
falls with time. Alternatively, if the drag term is negligible, 
then the clump will accelerate freely, and the velocity will be 
given by i = gt, if g is a constant. However, this case appears 
to be ruled out by observations showing slower velocities. 

Since material is continually stripped from the clump 
as it moves relative to the surroundings, it may eventu- 
ally disappear completely. Clearly, if this happens during 
a timescale that is less than that required for the velocity 
to exceed the terminal velocity, then we need not simulate 
the case of a clump moving at its terminal velocity. Assum- 
ing constant gravity, the free-fall time from 30 kpc is about 
10* yr for the value of g assumed above. If the clump started 
from that height with zero velocity with respect to the centre 
of the galaxy, it will have reached a speed of roughly 400 km 
s~^ in the free-fall time. Depending on the wind parameters 
and gravitational field, a 1O*M0 molecular clump will not 
reach terminal velocity and may not be totally evaporated 
away during free fall from 30 kpc. 

2.4 Origin of the clouds 

As mention above, the data suggest that much of the mate- 
rial within the filaments is outflowing. There are two scenar- 
ios in which filaments with the observed properties could be 
produced. Firstly, clouds may form in cooling gas accreting 
from outside the galaxy. As they fall inwardly they would be 
stripped by an outflowing wind, and the lower regions of the 
fllament could be infalling while the outer portion could be 
outflowing. Alternatively, such clouds may condense in gas 
that has been expelled from the central galaxy, and eventu- 
ally begin falli ng back into the centra l galaxy as in a galactic 
fountain (e.g. IShapiro fc Field! [l976l ). As the clouds fall in- 
wardly, they would be stripped by outflowing material. In 



2.5 Outflows 

In a previous subsection we discussed the idea of ram pres- 
sure due to an impinging wind, or outflow. However, it is 
unclear how precisely how fast these winds are: the density 
and emission structures of transonic flows can be similar to 
those of nearly hydrostatic media. So, in the absence of X- 
ray observations with high spectral resolution, direct obser- 
vational diagnosis of the kinematics is unable to distinguish 
flrmly between low Mach number and transonic speeds. De- 
spite this. X-ray maps do show evidence of weak shocks in 
the hot material at radii in the range of a few kpc to at 
least a few lO's of kpc in the Virgo and Pers eus clusters (e.g. 
iForman et al.ll2007l : ISanders fc Fabianll2007l '). Thus, we know 
that there are large-scale flows with transonic speeds. The 
sound speeds of 3 x lO^K and 10*K ionised gas, for the Virgo 
and Perseus clusters, are 8 x 10^ cms"'^ and 1.5 x 10* cms~^, 
respectively. In conjunction with the X-ray maps, this sug- 
gests that the wind speeds are likely to be ~ 10*cms~^. 
In fact, we suspect that the flows may be fountain- like and 
that outwardly flowing gas eventually falls back inwards. As 
a result, the wind velocity is not constant, but drops off with 
radius. So by the time it reaches the cold clouds the wind 
velocity will be much reduced from its central value. 

AGN activity is considered a prime candidate for pre- 
venting the ICM cooling from 10*K at rates of more 
than lOOMQyr"^. This requires an energy input rate of 
10'''*ergs~^ which is comparable to the kinetic energy flux 
of material with a mass density of 2 x 10~^® gem"'' flowing 
through a spherical surface of radius 10 kpc with a radial 
velocity of 10*cms~^. Thus, the "cooling problem" and our 
wind speed estimates imply comparable power inputs from 
AGN. 

The outflows we envisage would probably be somewhat, 
although not hugely, collimated. X-ray images of the centre 
of the Perseus cluster show some asymmetry in structures 
at radii approaching 10 kpc. Outside about 10 kpc, the in- 
teraction of AGN plasma with ambient galactic interstellar 
and intracluster gas appears to reduce asymmetries associ- 
ated with collimated AGN outflows. The filaments extend 
several times this far in radii, and are radial suggesting that 
the wind may be fairly isotropic at such radii. 

One would also expect the outflows to be somewhat 
time- dependent. However, the presence of clouds dotted 
throughout the ICM may well have the effect of obstruct- 
ing the outflow of energy injected by the AGN. Therefore, 
the time-variability of the AGN-driven flow far away from 
its source is likely to be fairly uncorrelated with the vari- 
ability at the source. Essentially this means that, although 
the initial energy injection rate may be time-variable, at 
larger radii the variations may have been smoothed-out mak- 
ing the flow relatively steady. Finally, intermittancy on long 
timescales much greater than 3 x lO'^ yr will be on timescales 
greater than the ages of the structures we are attempting to 
explain, and consequently unimportant in this work. Based 
on these considerations, we will assume the outflows to be 
transonic, isotropic and steady. 
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3 NUMERICAL SETUP 

3.1 Assumptions and methods 

We use hydrodynamical simulations to investigate the inter- 
action between a tenuous flow and a mass source. The source 
of mass is assumed to be a cold cloud from which material 
is stripped due to the interaction with an inci dent wind. 
The b a sic approach is similar to t hat adopted b y Falle et al] 
(|2002l ). IPittard et all (j2005l ). and iDvson et al] (|2006[ l. The 
main difference is that here we also include the effects of 
gravity to make the model more appropriate for a typical 
galaxy cluster environment. 

Despite including the gravitational potential, we will 
ignore self- gravity for this set of simulations, for the follow- 
ing reasons. Self-gravity is certainly important in the cold 
clump at the bottom of each filament, and would be impor- 
tant in our model if we were to try to calculate the mass loss 
rate of the clump. As we have already stated, the details of 
the injection process are not important for the question we 
are addressing. Self-gravity may well be important in the 
filaments, but the absence of a population of bright, young 
stars associated with the filaments (e.g. Hatch et al. 2006) 
also provides some evidence that it may not. As an illustra- 
tion, we give the ratio of the magnitude of the gravitational 
field to the magnitude of the thermal pressure force per unit 
volume in a cylindrical filament of mass density p and radius 
R is approximately, 



gpR 
P 



2nGpnmpR^ 



0.3 



R 



g cm~3 



30kpc 



10^ K 



(5) 



where prrip is the mean mass per particle in the gas. Thus, 
filament parameters are in the range that may allow the 
ratio to approach unity. However, it is clear that the values 
are sufficiently uncertain that an initial focus on cases in 
which self-gravity is neglected is reasonable. Furthermore, 
if self-gravity does lead to instability, shear may limit the 
development of clumps. We will neglect self-gravity for the 
present investigation. 

The simulations are split into two groups. In the first we 
assume constant gravity, and work in the frame of the cloud. 
This is to gain an understanding of the effect of the gravita- 
tional potential and the density contrast between the cloud 
and the wind. In the second group, we employ a realistic 
cluster atmosphere and work in the rest frame of the clus- 
ter, and advect the cloud. This allows us to implement what 
we have learnt from the prior simulations in the more real- 
istic case. The latter setup is discussed more fully in section 
6. In either case, the velocity of the incident wind increases, 
as measured in the clump frame, with time. Gravity also 
affects the motion of material in the tail. 

For all models we assume a rate of mass injection per 
unit time and volume that may vary temporally, but that is 
spatially uniform within a given radius of the cloud's centre. 
This mass injection prescription is simple to implement, but 
has the disadvantage that the flow from the injection region 
is isotropic. However, the effect of asy mmetric mass loss i s 
small at large distances from the cloud (jPittard et al.|[2005h . 
Therefore, we emphasize that the actual details of the mass 
injection process are relatively unimportant. The boundary 



of the mass injection region is not meant to coincide with the 
boundary of a cloud, and the cloud could be much smaller 
than the injection region. 

Three-dimensional calculations are required if the 
sources are spherical. To reduce the computational cost we 
restrict ourselves to two-dimensional simulations in which 
the sources are cylindrical. While we expect some differ- 
ences between calculations performed in two and three di- 
mensions, at this stage we can still gain important insight 
from less computationally demanding two-dimensional sim- 
ulations. 

In this study we treat the hot gas as adiabatic, since 
in the Perseus cluster the cooling time of X-ray emitting 
gas at 30 kpc is ~ 10^ yr, much longer than the flow time 
of the outflow along the tail. The cooling rate of the gas 
in the tail itself is very high, and must obviously be offset 
by a heating rate. We do not know the precise temperature 
of this gas, but the optically emitting gas in the tail is at 
about lO^K. An isothermal treatment of this material is rea- 
sonable. Thus, we investigate the simple case in which the 
incident wind behaves adiabatically, while the injected gas 
remains isothermal. 

To ensure isothermal behaviour, we use an advected 
scalar, a, which is unity in the injected gas and zero in the 
ambient gas. The source term in the energy equation is then 



K.ap{To - T), 



(6) 



where p and T are the local mass density and temperature, 
and where k is large enough that the temperature always 
remains close to the equilibrium temperature. To, in the in- 
jected gas. Inside the injection region we added an extra 
energy source so that th e gas is injected with temperature 
To (see lFalle et~alll200l for further details). 

The simulations are performed using mg.gt, a second- 
order accurate code with adaptive mesh reflnement (AMR) . 
mg^gt uses a hierarchy of grids G" — such that the mesh 
spacing on grid G" is Axo/2n. Grids G" and G^ cover the 
whole domain, but the flner grids only exist where they 
are needed. The solution at each position is calculated on 
all grids that exist there, and the difference between these 
solutions is used to control reflnement. In order to ensure 
Courant number matching at the boundaries between coarse 
and flne grids, the time-step on grid G" is Atg /2n where Ato 
is the time-step on G". 

In this model, we are studying the behaviour of cold 
clouds located ~ 10 kpc away from the centre of the cluster. 
At such distances, the gravitational acceleration in galaxy 
clusters is ~ 10~^cms~^. We investigate different density 
contrasts between the cloud material and the ambient gas, 
and compare the morphology of the results with observa- 
tions. 



3.2 Variables 

The problem we are investigating is scale-free. For numer- 
ical reasons, we have chosen scales so that the values of 
parameters are near unity. For easier comparison with the 
astrophysical systems of interest we will rescale the results 
to physical units. In particular, we have scaled the length by 
a factor L, the gravitational acceleration by G, and the den- 
sity by R. In this case, the hydrodynamic quantities scale 
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X — Lx' , 

9 = Gg, 
p = Rp, 

t = Tt' 



L , 

V = —V , 

r 



P ■ 



R 



P ■ 



(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



where the primed quantities represent the dimensionless 
simulation values, and the unprimed quantities represent the 
real values, t is time, v is velocity, and p is pressure. In cases 
in which q is nonzero 



(13) 



In our simulations the values selected for parameters are 
roughly consistent with estimates of the properties of clouds 
in the Perseus cluster: the length L = 50 pc, G = 10~^cms~^ 
(except for one case) and R = 10~^^gcm~^. For cases in 
which <; = 0, we leave r — 1.24 Myrs, which is the same 
value as when g is nonzero. In section 5, we have employed 
more accurate values of the scaling factors. 



3.3 Initial conditions 

To ensure that the ambient wind does not penetrate the 
injection region, we take the ram pressure of the injected 
material at the boundary of the injection region to be larger 
than that of the wind, i.e. we require. 



2 

Vi pi 



P^Vv. 



(14) 



where ri is the radius of the injection region, Vi is the flow 
speed of the injected material at r = ri and pi its density 
at this point. Q is the rate per unit volume per unit time at 
which mass is injected within r — Vi. 

If we use the pressure behind a stationary normal shock 
in the wind rather than the wind's ram pressure, the above 
criterion becomes 



Pv/V^fi 



7+1 ViTi 



(15) 



In the frame of a falling clump, the wind's ram pressure 
increases with time. Thus, the minimum mass injection rate 
must increase accordingly to ensure that the wind does not 
penetrate the injection region. If the value of Q at t = is 
Qo and the gravity is constant, the ratio of Q to the wind 
ram pressure will remain constant if 



1 + 



9t 



(16) 



where Ww,o is the relative velocity between the clump and the 
ambient flow at f = 0. While the method that we are using 
to treat mass injection requires that we take Q to increase 



with time, in reality t he mass injection rate w ill increase 
with ram pressure (e.g. iHartauist fc Dvsonlll988l ). 

The computational domain is —200 ^ a; 200, —100 5S 
y ^ 400, with the mass injection region centred at the origin. 
We use 6 grids, and the size of the finest cells is 0.125. 

By scaling the physical variables we can choose units 
such that ri = 1, = 1, and p^ — 1 and initially set = 
10 and To — 1. Furthermore, we assume that the injected 
material is in pressure equilibrium with the ambient wind. 
Thus, if we assume that the injected material is a factor of 
10'^ cooler than the ambient wind, then it is also 10'^ times 
denser than the wind. We will use rj to represent this contrast 
parameter. 

We investigate the morphology of tails for three differ- 
ent contrast parameters: rj — 10^, 10'^, and 10*. The corre- 
sponding initial values of the adiabatic Mach number, M, 
are: 0.77, 0.24, and 0.077, respectively. These parameters 
define case 1), 2), and 3), respectively. Clearly, as the am- 
bient wind accelerates with respect to the cloud, the Mach 
number will increase. 

For comparison, we also study several other cases in 
which the density contrast is 10'^. Case 2i) is identical to case 
2) except that g' — 10 and G — 10~*. In cases 2ii) and 2iii) 
the gravity is zero, but there is a constant wind velocity of 10 
and 30 (in code units) respectively, corresponding to Mach 
numbers of 0.24 and 0.72. These parameters are summarised 
in table 1. It does not make sense to specify the mass of 
the clouds, since they are likely to be composed partly of 
molecular material, and therefore considerably more massive 
than an estimate of 47rrj?pc/3 would suggest. 

The code wind velocity v'^ — 10, corresponds to a ve- 
locity of iiw = 3.9 X lO^cms"^, which, is comparable with 
the transonic values described in the previous section. 



3.4 Sub-grid turbulence model 

Simulations with a transonic wind pose difficulties that are 
absent in those for hypersonic winds. Firstly, in order for the 
boundaries to have no effect on a simulation with a transonic 
wind, the computational domain has to be very large. Sec- 
ondly, the velocity shear between the injected material and 
the wind is so extreme in the transonic case that the wind 
flow separates and produces a turbulent wake downstream 
of the interaction region. Calculations based on the Euler 
equations cannot adequately describe such turbulence, and 

a turbulence model must be employed. 

We adopt the k — e model used by iFalle et al.l l|2002l ). 
The subgrid turbulence model allows the simulation of a 
high Reynolds number flow through the inclusion of equa- 
tions for the turbulent energy density and dissipation rate, 
from which viscous and diffusive terms are calculated. The 
model solution gives an approximation to the mean flow. 
The diffusive terms in the fluid equations model the turbu- 
lent mixing of the injected gas with the original flow due to 
shear instabilities. The model has been calibrated by com- 
paring the computed growth of shear layers with experi- 
ments (Dash & Wolf 1983). It is based on the assumptions 
that the real Reynolds number is very large, which is the 
case in astrophysical flows, and that the turbulence is fully 
developed. The use of the subgrid model gives more realis- 
tic results than those obtained from a hydrodynamical code 
where the sizes of the shear instabilities are determined by 
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Table 1. Simulation parameters: ambient density, cloud density, sound speed of injected material, initial 
Mach number of the wind velocity with respect to cold material, gravitational acceleration. The implied 
mass outflow rates close to the origin are i-wr^ 5p^^^|^v^ ^ lOM0yr— where <5 ^ 0.1 is the fraction of the 
sky covered by the outflow. 



Case 


Pamb (9 cm 3) 


pcigcm 3) 


Cs,c(cms ^) 


Mnitial 


9 (cms 2) 


1 


10-26 


10-24 


5.2 X 10^ 


0.77 


10-7 


2 


10-26 


10-23 


1.6 X 10^ 


0.24 


10-7 


3 


10-26 


10-22 


5.2 X 10* 


0.077 


10-7 


21 


10-26 


10-23 
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the numerical resolution. Further details can be found in 



iFalld l| 19941 ). 



4 RESULTS 
4.1 Morphology 

In general, the morphologies of the filament density gener- 
at ed by accele rating fiows are very similar to those presented 
in lPittard et a l. (2005). A bow wave develops upstream of 
a mass source, and the wave's amplitude falls off as 1/r. 
If the flow becomes supersonic, the bow wave becomes a 
bow shock and a very weak tail shock may occur in the 
wind downstream of the mass source, when the external flow 
is subsonic. The injected material remains in approximate 
pressure equilibrium with the wind, and eventually becomes 
confined to a tail with a width comparable to that of the 
injection region (cf. figures [1] to |4]). The opening angle of the 
injected material increases with the external Mach number. 

It should be noted that we use a sub-grid turbulence 
model which, in effect, models a turbulent boundary layer 
at the surface of the filament. So even when the tail appears 
smooth there is significant turbulence. For increasingly large 
density contrasts, the turbulent viscosity, which is an inte- 
gral part of the subgrid turbulence model, fails to stabilize 
the subsonic shear layer. In this case, there is turbulence on 
the grid as well as sub-grid turbulence. Large-scale insta- 
bilities in thin filaments demand high density contrasts and 
an ambient flow that remains subsonic with respect to the 
clump. 

The opening angle of a turbulent boundary layer is 
given by Canto & Raga (1991) as. 



dh 
dx 



2 



(17) 



where a is the so-called spreading parameter. Experiments 
show that at low Mach number, turbulent mixing layers have 
a constant opening angle, and a ~ 11, so that dh/dx ~ 10.4. 
At high Mach numbers, a ~ 50, so dh/dx ~ 2.3. The tail 
becomes fully turbulent at a downstream distance of. 



CO 



2e 



(18) 



where e is the entrainment efficiency, Mw is the Mach num- 
ber of the wind with respect to the clump, Ca,a is the sound 
speed of the environment, Cs,c is the sound speed of the tail 
material, rc is the initial radius of the tail. 



To calculate this length. Canto 
16.9A'fwrc, which can be rescaled as. 



Raga use x 



X ~ 0.85 kpc 



1 



50 pc 



(19) 



in other words, for sensible values of the wind fiowing past 
a clump, we should expect the tail to be turbulent on kpc 
scales. 

These figures show a clear dependence, at least at early 
times, of the tail morphology on the density contrast be- 
tween the cloud and the ambient gas. For example, for 
rj = 10* (figure 3), the tail is nearly linear, but is punc- 
tuated by fairly evenly spaced blobs along the length of the 
tail. This is because subsonic shear layers tend to be unsta- 
ble unless the viscosity is high. Some turbulent jets display 
similar behaviour. The blobs are separated by increasingly 
large distances further along the tail, indicating accelera- 
tion. For 77 = 102, (figure 1) the tail only becomes irregular 
at large distances away from the cloud, while for rj — 10'^ 
(figure 2) most of the tail is irregular. For comparison we 
also show the morphologies for a case when the gravity is 
weaker (figure 4) and nongravitational cases with different 
ambient flow velocities (figures [5] and [6]). 

For ambient fiows that remain subsonic in the frame 
of the clumps, the tails are irregular in appearance. In case 
2i) the gravity is strong enough that the initially subsonic 
flow becomes supersonic, and consequently the tail becomes 
smooth and regular. Morphological considerations by them- 
selves would suggest that the filaments in the Perseus clus- 
ter satisfy these conditions. However, an examination of the 
kinematic structures of the model and observed filaments is 
also necessary. 

4.2 Velocity profiles 

Figures [7] to [10] show the mean velocity of the tail material 
at 100 increments along the length of the tail for each of 
cases 1), 2), 3), and 2i). 

Except for case 2i) and to some extent case 1), each 
velocity profile shows some scatter around a mean velocity 
that is roughly constant along the length of most of the fil- 
ament. We do not show velocity profiles for cases 2ii) and 
2iii) because they show behaviours similar to those for cases 
2) and 3). However, the velocity profiles, including those 
for the model filaments that bear the most morphological 
resemblence to the observed filaments, are not compatible 
with the observed velocity structures of the NGC 1275 fila- 
ments. Nevertheless, the third panel of figure [10] shows clear 
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Figure 7. Velocity profile for case 1) for (from left to right) t' = 6.95, 14.5, 20.9, and 28.2 (t = 8.3, 18, 26, and 35 Myrs), respectively 
The Mach numbers of the ambient flow with respect to the cold clump are: 1.3, 1.9, 2.4, and 3, respectively. We have removed the central 
10 kms~^ to avoid any points where the tail does not exist. 
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Figure 8. Velocity profile for case 2) for t' = 3.4, 8.3, 12.1, and 16.9 {t = 4.3, 10, 15, and 21 Myrs) respectively. The Mach numbers of 
the ambient flow with respect to the cold clump are 0.33, 0.45, 0.54 and 0.66, respectively. We have removed the central 10 kms^^ to 
avoid any points where the tail does not exist. 




C :i . C J j . 3 




Figure 1. Density morphology for case 1) at t' RJ 28 (t ~ 35 
Myrs) at which time the external flow at the position of the clump 
has accelerated to M = 3. The initial Mach number was 0.77. Each 
spatial unit corresponds to 50 pc. Thus, the plotted domain, in 
the x-direction, extends between ± 1.5 kpc. The density contrast 
between the cold clump and the ambient gas is lO'^. Note the bow 
shock around the clump, and also the tail shock behind this. 



zig-zag features, which are similar to t hose in the v elocit y 
profiles given in figures 5, 7 and 14 of iHatch et al.l (|2006l ). 
Note that the zig-zag features in figures [7] and [10] also corre- 
spond to the regions of enhanced density along the tail. The 
similarity between figures [71 and llOl suggests that the zig-zag 



Figure 2. Density morphology for case 2) at t' 17 (t ~ 21 
Myrs) when M = 0.66. The initial Mach number was 0.24. The 
density contrast between the cold clump and the ambient gas is 
10'^ . Note the irregular structure of the tail along the entire length 
of the filament, compared to the linear tail in the previous figure. 



features occur when the Mach number of the ambient flow, 
with respect to the clump, is roughly 2-3, after ~ few 10s 
Myrs. 

The velocities near the heads of the model filaments are 
low, whereas those of the observed filaments are not. The 
discrepency could be due to the material near the heads of 
the observed filaments being warm, but neutral, and hence 
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Figure 9. Velocity profile for case 3) for t' = 1.93, 5.7, 12.7 and 18.1 (t =2.4, 7.0, 12, and 22 Myrs) respectively. The Mach numbers of 
the ambient flow with respect to the cold clump are 0.09, 0.12, 0.15, 0.18 and 0.22, respectively. We have removed the central 10 kms~^ 
to avoid any points where the tail does not exist. 





Figure 3. Density morphology for case 3) at t' 18 (t ~ 22 
Myrs) when M = 0.22. The initial Mach number was 0.077. The 
density contrast between the cold clump and the ambient gas is 
10''. Note the approximately linear tail, with at least 6 descrete 
blobs of material close to the clump. 



Figure 4. Density morphology for case 2i) at t' 20 (t ~ 25 
Myrs) when M = 5.2. The initial Mach number is 0.24. The den- 
sity contrast between the cold clump and the ambient gas is 10^. 
The structure of the tail in this example is regular, with no obvi- 
ous blobs at various locations along the tail. 



invisible in the observed optical lines. More likely, it is a re- 
sult of our prescription for mass injection which will produce 
different features near the injection region, but the large- 
scale downstream flow should be unaffected. 

The fact that the profile shown in figure [TO] resembles 
the observed profile better than that shown in figure [71 
is probably because the density contrast between the cold 
clump and the ambient flow is more realistic in case 2i) . It is 
likely that if we had continued simulating cases 2) and 3) to 
high enough Mach numbers that they may also have begun 
to exhibit the zig-zag structure. Based on these results it 
appears that the velocity fluctuations (and scale-size of in- 
stabilities) are larger for greater density contrasts between 
the clump and the wind. 



5 TAILS IN A REALISTIC CLUSTER 

ENVIRONMENT: THE PERSEUS CLUSTER 

Having investigated the properties of accelerating flows past 
cold clouds, we now study the more specific case of a galaxy 
cluster atmosphere. Firstly, we assume a generic gravita- 



1 



Figure 5. Density morphology for case 2ii) at t' ft: 20 (t ~ 25 
Myrs). The density contrast between the cold clump and the am- 
bient gas is 10^, and gravity is omitted from the simulation. M 
= 0.24. 
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Figure 10. Velocity profile for case 2i) at t' = 2, 5, 10, and 20 (t = 2.5, 6.2, 12, and 25 Myrs), respectively. The Mach numbers of the 
ambient flow with respect to the cold clump are 0.74, 1.5, 2.7 and 5.2, respectively. We have removed the central 10 kms~^ to avoid any 
points where the tail does not exist. 




Figure 6. Density morphology for case 2iii) at a code time of 
t' 20 (t ~ 25 Myrs). The density contrast between the cold 
clump and the ambient gas is 10''. Gravity is omitted from this 
simulation. M = 0.72. 



tional potential which leads to a /3-profile gas density profil e 
in hydrostatic equilibrium (e.g. IVernaleo fc Reynold3l2006l ). 



P = 



Po 



(20) 



where we have assumed spherical symmetry, r is the radial 
coordinate, and ro is the scale-height of the density distri- 
bution. The gravitational acceleration is therefore given by 



(21) 



where T is the gas temperature which is assumed to be con- 
stant with radius, ^ — 0.6 is the mean mass per particle. 

For each such model, we will perform the calculations 
in the frame of the central galaxy and prescribe how the 
cloud injection region moves. Due to the structure of the 
gravitational potential employed and the initial location as- 
sumed for the cloud, the solution of the cloud motion is very 
close to that for free-fall. Assuming free-fall, we find that the 
velocity, after starting from rest at xo, is. 



v{r) — - 






'■f) + Xq 


1/2 


(22) 















Equation (|22|) can be numerically integrated, using the 



values for specific galaxy clusters, to obtain the position of 
the cloud at each point in time, r(t). For example, for the 
Perseus cluster we use, T = 10* K, po = 7 x 



lO^^^gcm-3 



P = 0.81, and ro = 28.5 kpc (e.g. iPope et al.ll2006l ). With 
these values, the motion of a cloud falling freely from 30 kpc 
can be well approximated by, 



r{t) = 30 



0.9 



t 



75 Myr 



kpc. 



(23) 



which is roughly what one would expect for constant gravity. 
This correspondence occurs because the gravitational accel- 
eration scales as 1/r at large radii, and r at small radii. The 
cloud is initially located around the turning point where the 
function is relatively flat. This means that constant gravity 
is a suitable approximation and the previous set of simula- 
tions, in which constant gravity was assumed, retain their 
validity. 

For completeness, the Mach number of the flow around 
the clump, with respect to the ambient flow, is 



M 



2/3 



7 



1/2 



(24) 



Thus, for /3 — 3/2, ro = 28.5, xq = 30 kpc, we flnd that 
M > 1 for r < 0.9ro ~ 26 kpc. Consequently, for higher 
wind velocities, the flow becomes supersonic closer to the 
initial location, xq. The ram pressure due to relative motion 
between the cloud and the static atmosphere is. 



p[r{t)] 



dr{t) 



dt 



(25) 



We choose an initial density contrast between the cloud 
and its surrounding of 10*, and set the scaling parame- 
ters to i? = 7 X IQ-^^gcm-^* and L = 3.085 x 10^^ cm. 
P = pk\^T/fj,mp « 3.6 X 10"^" erg cm""' in the centre 
of the Perseus cluster, giving G = P/{P'RL) = 1.66 x 
10~^°cms~^. The temporal scaling constant, is r = 4.3 x 
10^^ s, while the velocities are scaled by L/r — 7.2 x 
10® cms"'^. The time taken for the cloud to fall from 30 kpc 
to 3 kpc is 75 Myr, or 0.55 in terms of code time units. For a 
density contrast of 10^, we use R = 7x 10~^^ g cm""' and L — 
3.085 X 10^^ cm, giving G = 1.66 x 10"^ cms'^ This means 
that the temporal scaling constant becomes r = 1.4 x 10^® s, 
and the velocities are scaled by L/r = 2.2 x 10® cms~^. The 
time taken for the cloud to fall from 30 kpc to 3 kpc now, 
or 1.7 in terms of code time units. 

Finally, in the preliminary simulations we found that 
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the tail widths were comparable to the diameters of the 
clouds. Therefore, to generate filaments of the observed 
widths which are constrained to be <1 kpc we will simulate 
the effect of material stripped from relatively large clouds, 
with an injection region of radius 100 pc. As already stated 
this may be significantly larger than the physical size of the 
cloud. 



5.1 Results 

Since the lateral velocity dispersion (zig-zag effect) appears 
to increase for larger density contrasts, we will concentrate 
on the case r} = 10*. However, the case where rj — 10^ is also 
of interest. The general properties are as follows: long tail- 
like structures only form if the wind velocity is sufficiently 
large. This wind is assumed to be the same material, at the 
same temperature as the static atmosphere. For example, for 
low values of v'^ (5 and 10, corresponding to 3.6 x 10^ and 

7.2 X 10 cms^ for 77 = lO"*), the tail never becomes long 
because the ram pressure due to the impinging wind cannot 
overcome the gravitational force acting on the tail material. 
As a result any tail that forms is infalling everywhere, which 
contradicts the observational data, at least for NGC 1275. 
It therefore appears that the length and morphology of the 
tails are extremely dependent on the wind velocity. 

The morphology of the tail also depends strongly on the 
density contrast between the cloud and its environment. In 
each figure, the cloud is shown as red and is the lowest point 
along the line defined by x=0. Material is injected (stripped) 
and collects behind the free-falling cloud, and some dense 
material sinks either side of the cloud. Behind this cold ma- 
terial, a filamentary tail forms. Inhomogeneities in the tails 
are dependent on the wind speed and density contrast, in 
agreement with the previous results. As can be seen from 
figures [TT1 to ll4l it appears that a density contrast of 77 = 10^ 
produces structures most qualitatively similar to those seen 
in real filaments. However, they are low density, being a few 
X 0.1 particles per cubic centimetre. This does not present 
problems in terms of pressure balance between the tail and 
the ambient medium since the tail material is diluted and 
heated to temperatures in the range 10^ - lO'^K. The simula- 
tions produce filaments of 10*' — 10^ M0 while the estimates 
in the following section suggest that 10* M© filaments might 
be more likely. Even so, the quantity of cold gas, in the sim- 
ulated tail, is more th an capable of prod ucing the observed 
optical luminosity (see lHatch et al.ll2007l . for example). 

The material around the head of the cloud also exhibits 
different behaviour in this set of simulations. Material is 
stripped from the cloud by the wind and builds up behind. 
It appears that in many cases much of this material does 
not make it into a tail but simply collects behind the cloud. 
Eventually this becomes unstable, probably due to Rayleigh- 
Taylor instabilities, and produces the structures seen in the 
simulations. Similar structures have b een noticed in previou s 
numerical simulations in the work bv iMurrgiv fc LinI (|2004 ). 

In the simulations, initial wind velocities of ~ 
10* cms~^ are required to produce long tails which are out- 
flowing at their outer extent. This large velocity is expected, 
since the velocity of the wind, at the clump, must exceed the 
infall velocity of the clump for material to be carried away 
from it. By the time the wind has reached the clump it has 




Figure 11. The density structure at 58 Myrs, for r) = 10 and 
wind velocity = 5.8 X lO^cms"^. Panel shows data at 3.4 X 
10^ yrs. The length scale indicates the distance from the cluster 
centre in kpc. The density scale shows the density in code units 



''1 ' 



Figure 12. The density structure at 38 Myrs, for rj = 10* and 
wind velocity = 4.32 X 10* cms~^. The length scale indicates the 
distance from the cluster centre in kpc. The density scale shows 
the density in code units (p'). 



done considerable work against gravity. This leads to the 
velocity dispersion along the clump of a few xlOOkms"^. 

Figure [T5] shows the velocity profile for 77 — 10* and a 
wind speed of 5.8 x lO^cms"^ (i;4, ~ 80). This case proba- 
bly best resembles the observations. Figure[TS]seems to show 
that a lower wind velocity (u^, = 60, Uw = 4.32 x 10^ cms~^) 
does not produce the observed zig-zag kinematics for this 
particular density contrast. These results should be com- 



pared with the 77 = 10 case. 

Figures [17] and [18] show the kinematics for two wind 
velocities: u(^, — 80 and 60, in figures [17] and [18] These cor- 
respond to 1.76 X 10*cms"^ and 1.32 x 10*cms"^ which 
are either side of the sound speed of the ambient medium 
(1.52 X 10*cms~^). Figure [18] eventually shows kinematic 
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Figure 15. Velocity profile for case rj = 10'* and wind velocity = 5.8 X lO'^cms *. Panels show data at 14, 21 and 34 Myrs. The length 
scale indicates the distance from the cluster centre. We have cut off the central 20 kms^^ to exclude points where there is no tail. 




Figure 13. The density structure at 20 Myrs, for r} = 10'^ and 
wind velocity = 1.76 X 10* cms~*. The length scale indicates the 
distance from the cluster centre in kpc. The density scale shows 
the density in code units (p'). 



Figure 14. The density structure at 24 Myrs, for rj = 10^ and 
wind velocity = 1.32 X 10* cms~*. The length scale indicates the 
distance from the cluster centre in kpc. The density scale shows 
the density in code units (p'). 



structure that could be interpreted as a zig-zag, although it 
is not as clear as in figure 1151 The tails extend to greater 
lengths for rj — 10"^ than rj = 10* because the evaporated 
material is less dense and can be more easily transported 
outwards by the wind. 

Again, the amplitude of the zig-zags is smaller for 
rj = than 77 = 10*. As a result, these simulations suggest 
that a density constrast between the cold material and the 
ICM of 10* is a good description, and that a wind of ve- 
locity ~ 6 X 10^ cms~* is required to produce the observed 
filamentary structures. Furthermore, the kinematic zig-zag 
features only seem to form after about 30 Myrs, suggesting 
that the Perseus filaments are approximately this old. 

In each case, the total mass removed from the cold cloud 
was of the order of 10^ Mq . However, the majority of this 
material is not in the filament, but is in the more amorphous 
and disrupted region behind the cloud. The material around 
the cloud is likely to have a low ionisation and is therefore 
not likely to be optically emitting. 



6 MOMENTUM TRANSFER FROM THE ICM 

The kinematics of the filaments indicate that momentum 
transfer between the cool and hot phases of the ICM occurs. 
If there is a large enough quantity of cool gas, this could 
have important consequences for the flow of the hot phase. 
Momentum P is transferred from the hot ICM phase to a 
cloud at a rate 

Pciump = KAp„{v^ - x)'^, (26) 

If the radius of the cloud is ~ 30pc, and (uw — i) ~ 
10'^-*cms-* then Pciump ~ lO^^-^^ergcm"*. 

A rough estimate for the rate at which momentum is 
transferred from the external medium to the cloud's tail is 

-Ptail = "Itailitail, (27) 

where mtaii = Apt^ultau , Ptaii is the density in the tail, itaii 
is the typical acceleration in the tail, and Ztaii = k is the 
length of the tail. Constant cross-sectional area, density, and 
acceleration along the tail are assumed. The acceleration can 
be approximated by itaii = (i — ii)^/2Ztaii where x is the 
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Figure 16. Velocity profile for case rj = 10'* and wind velocity = 4.32 X 10** cm s ^ . Panels show data at 16, 23 and 38 Myrs. The length 
scale indicates the distance from the cluster centre. We have cut off the central 20 kms^^ to exclude points where there is no tail. 
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Figure 17. Velocity profile for case rj = 10^ and wind velocity = 1.76 X 10** cms ^. Panels show data at 8, 12 and 20 Myrs. The length 
scale indicates the distance from the cluster centre. We have cut off the central 20 kms~^ to exclude points where there is no tail. 



velocity of the clump, and u is the maximum velocity of the 
tail material. If the number density in the tail is about 10 
cm~^, the tail radius is about 50 pc, and the magnitude of 
\u — ij is ~ 10^ km s~^, then 

^*^"~'°''(l0^)(5fe) ( lOo\m"s-0 

(28) 

The corresponding rate at which systemic kinetic energy is 
transferred to the tail is iStaii ~ -PtaiiK^t ^ 

[w^)[50Vc) l lOOkms-i J 

(29) 

If the energy dissipated as thermal energy were comparable 
to this and all of that thermal energy were radiated in the 
observed optical lines, about 100 tails would account for the 
observed optical luminosity of the NGC 1275 filaments. 

The estimate above suggests that the tails are extremely 
massive ~ 10* M0 , t hough the mass of i onised material must 
be much lower (e.g. iHatch et al.ll200^ . This suggests that 
the filaments may consist of molecular material surrounded 
by a skin of ionised material. 

We can compare the momentum transfer between the 
phases of the ICM, with that developed by an AGN jet, 
which may in turn drive a wind from the central galaxy. A 
jet injecting material at lM0yr~^ with a velocity of 0.1 



times the speed of light over a lifetime of 10 Myrs, has 
P ~ 10^°gcms~^. This can be compared to the total mo- 
mentum transferred along the length of a filament by the 
wind ~ iAitaiiptaii ~ 10**gcms~^. This suggests that the 
momentum transfer to many tails could be important for 
dissipating the energy injected by an AGN in the ICM. In 
particular, since much of the cold material is located near the 
centre of the cluster, if this dissipation mechanism is signifi- 
cant, then it will increase the AGN heating efficiency in the 
central regions of the cluster. Without this obstructing cold 
material, more of the injected energy would be deposited at 
larger radii, thus failing to heat the central regions where 
energy is most needed. 

To summarise this point it is worth reiterating that tails 
only seem to form if there is a significant outfiowing wind 
- they are formed by the ambient wind dragging and accel- 
erating material, transferring energy and momentum to the 
cold material. As a result they are probably a sign of activity 
in the central galaxy. This mechanism can provide a way of 
coupling the wind to the cold material. Given that the cold 
clouds are likely to proliferate near the cluster this means 
that the energy will be dissipated more centrally than in 
the absence of the cold material. As such, the cold material 
increases the effective opacity of the ICM to the outflowing 
wind. 

To estimate this effect, let us assume that energy is 
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Figure 18. Velocity profile for case r] = 10^ and wind velocity = 1.32 X 10** cm s ^. Panels show data at 0.9, 1.4 and 2.4 X 10^ yrs. The 
length scale indicates the distance from the cluster centre. We have cut off the central 20 kms^^ to exclude points where there is no tail. 



injected at a rate Ein, the energy dissipation rate per tail 
is i?taii, and there are A'^ tails. The change in power of the 
wind, over a region in which there are A*' clouds, will be 
AE = Ein — NEtaii- Taking the derivative with respect to 
radius, and assuming that -Etaii is a constant we find. 



dE _ dN ■ 

— — — -C'tail, 

dr dr 



(30) 



where dN/dr is obtained from the spatial distribution of the 
clouds. 



— — = 47rr n(r) 
dr 



(31) 



and n(r) is the concentration of clouds per unit volume. If 
we assume that the clouds are distributed similarly to the 
hot gas, in a /3-profile, then 

no 



(32) 



(l + (r/ro)2)/'' 

where no is the central density and ro is a scale length. 
Choosing [3 = 1 means that equation (|3Up can be integrated 
analytically to give. 



E{r) = Ein - NoEt^ 



r 
ro 



arctan 



(33) 



where A'o = 4nnoro- Therefore, the injected energy would 
be dissipated over a typical length-scale 

rdiss ~ ■■ • (34) 

Given that ro is probably a constant of the system, and 
Stall is also probably fairly constant, the important pa- 
rameters are -Ein and A^o- For a moderate AGN outburst, 
Ein ~ 10''^ ergs"\ with Staii ~ 10""' ergs"^ and A'o ~ 1000 
(corresponding to a total cold mass of 10^^ Mq), the dissi- 
pation length is of the same order as the scale-height of the 
cloud distribution. (This may be comparable to the scale- 
height of the ICM, but that is not clear. It is also possible 
that the clouds are more centrally concentrated than the 
hot ICM, but again, this is speculation.) For more powerful 
outbursts, a significant quantity of energy may be trapped, 
but the majority will still escape outside ro. 

It is worth noting that if there has been significant re- 
cent cooling (i.e. a deficit of heating) we may expect A'o to 
be larger, for a given system. For a given energy injection 



rate, the energy will necessarily be dissipated closer to the 
centre than in an otherwise identical system where A^o is 
lower. As a result, the presence of cold gas may provide an 
additional feedback process that helps to trap energy in the 
central regions. This effect will be greater if there has been 
a recent deficit of heating, with the trapping effect being 
reduced in the absence of cold material. 

This short derivation is, of course, a drastic simplifica- 
tion of real life. For example, in reality, the power dissipated 
in a tail is a reflection of the power of the incident wind, 
which will undoubtedly be a function of radius. However, 
it highlights the possibility that the production of filaments 
may be able to trap a fraction the energy injected by a cen- 
tral AGN. 



7 GENERAL EXPECTATIONS FOR OPTICAL 
EMISSION-LINE FILAMENTS IN GALAXY 
CLUSTERS 

The current study is centred on the Perseus cluster and con- 
cerns the interaction of material injected from one cloud into 
a hot flow unaffected by other clouds. A global consideration 
of the effects of a distribution of clouds on the hot flow out 
to ~10 kpc would be required to investigate the problem 
more completely, and we plan to investigate such problems. 
However, a full answer would require the developement of a 
detailed theory for the spatial and mass distributions of the 
clouds; given the state of the theory of our own interstellar 
medium this task would be very difficult. 

The nature and number of filaments will depend on: a) 
the spatial, mass and velocity distributions of cold clouds at 
the onset of AGN activity; b) the rate at which the AGN 
deposits mechanical and thermal energy into its outfiow and 
the asymmetry of these AGN outflows; c) the pressure, ve- 
locity, and density distributions of the hot intracluster gas 
at the onset of AGN activity; d) the gravitational potential 
of the galaxy. Below we consider, in sequence, how each of 
these could affect the fllamentary structure. 

a) If all other factors were the same, the total number 
of fllaments and their combined optical luminosity would 
be roughly proportional to the number of clouds that could 
survive interactions with an outflow long enough for flla- 
ments to evolve. If the distribution of cloud masses were 



14 E.C.D. Pope, J.M. Pittard, T.W. Hartqmst, S.A.E.G. Falls 



identical in all clusters, the systems with the most cold gas 
would exhibit the most luminous optical emis sion fr o m fila - 
ments. This is consistent with the findings of lEded l|200lh . 
The clouds which survive for longest are likely to be the 
most massive. 

b) The number of massive clouds will depend on AGN 
activity. If the clouds are formed in a steady fountain flow, 
the rate at which material turns into clouds would be pro- 
portional to the outflow rate of hot material as measured 
near its the source. Roughly speaking, this mass outflow 
rate may be proportional to the mechanical luminosity of the 
AGN wind. In addition, the destruction rate of the clouds 
may also be proportional to the outflow rate. For example, 
a higher outflow rate could lead to a more rapid destruction 
rate of the clouds. Thus, it is possible that the number of 
clouds producing fllaments would be fairly insensitive to the 
AGN mechanical luminosity. However, it would be cavalier 
to assume the locations at which clouds form, and the ini- 
tial formation of the clouds, to be independent of the outflow 
rate. 

c) The properties of the ICM could affect the initial 
mass function of the clouds. If signiflcant, thermal conduc- 
tion might establish the minimum lengthscale on which ther- 
mal instability grows. In hotter clusters, the minimum cloud 
mass would be higher, leading to a greater fraction of the 
clouds being formed with high masses. This would tend to 
increase the number of long fllaments. 

High ICM densities will lead to more rapid cooling and 
possibly a higher formation rate of clouds, and hence more 
fllaments. However, given that AGN feedback is likely to 
be important, high cooling rates may cause strong AGN 
activity and thus reduce the cloud formation rate. The exact 
balance between heating and cooling seems to vary from 
cluster to cluster. 

Clouds in the ICM may be so plentiful in the centre of 
a cluster that a bubble blown by the AGN becomes trapped 
near the cluster centre, preventing the formation of long 
fllament. So, high ICM density need not necessarily result 
in the production of many filamentary tails. 

d) The gravity of the central galaxy may play a role in 
limiting the size of the bubble mentioned in c). For example, 
the greater buoyancy may cause a bubble to detach sooner. 
Strong gravity may prevent an extended fountain or region 
of outflows. Even if it does not, it may limit the lifetimes of 
clouds as it would cause them to fall more rapidly. 

In summary, the question about scaling concerns very 
complex issues. A very long concentrated program would be 
required to address them properly, and it is unclear how 
reliable the results of such a program would be. 



8 SUMMARY 

This work is intended mainly as a preliminary study into the 
processes that shape fllamentary structures in the ICM, and 
not an exhaustive study. It is likely that there are many other 
processes such as magnetic fields, thermal conduction and 
viscosity that all play a role in the dynamics that produce 
the fllamentary kinematics and morphology. However, we 
have concentrated on the basic hydrodynamic system and 
have found several interesting mechanisms that may go some 
way to explaining the properties of the observed fllaments. 



The flrst ensemble of simulations show that the long, 
relatively straight fllaments at the centre of the Perseus clus- 
ter may be formed by the interaction of a wind with clumps 
of cold material. Flows of low Mach number winds generate 
structure in the filaments which may be comparable with 
observations. These results also suggest that the amplitude 
of the velocity fluctuations along the length of the filament 
grow with increasing density contrast between the cold ma- 
terial and the ambient wind. The morphology of the ob- 
served filaments suggests that the density contrast is large 
(~ 10*). Interestingly, the fiuctuations evolve with time pro- 
viding a possible diagnostic tool for estimating the ages of 
fllaments from observations, in this set of simulations, the 
best comparison with observations suggests fllament ages of 
~ 40 Myrs in the Perseus cluster. 

In a more realistic environment with spatially varying 
gravity, based on the Perseus cluster, fllaments only form if 
the wind velocity is sufficiently large. If there is no significant 
wind, then the optical emission may be more amorphous, or 
more likely, fllamentary on much shorter length scales. The 
morphology and kinematics suggests that a density contrast 
of 77 ~ 10^ niay be more compatible with the observations, 
while the best agreement in the kinematic data occurs for 
■q ~ 10*. As a compromise it is possible that the real condi- 
tions lay somewhere between these two extremes. Alterna- 
tively, other physical processes could be important. It should 
also be noted that these were 2-d simulations, and a full 3-d 
study would really be required to flnd the model conditions 
that best match reality. Regardless of this, the results do 
point to a relatively narrow range of parameters providing 
possible information about the state of the ICM, and the age 
of the fllaments. These results suggest that the cold material 
is of the order of 10* times denser than the ambient ICM, 
there is an outflowing wind of velocity of ~ 10^ kms~^, and 
that the fllaments are approximately a few lO's of Myrs old. 

The masses, and densities, of the fllaments created by 
these simulations seem to be low compared with values pre- 
sented in the previous section. The simulations produce fil- 
aments of 10^ — 10^ M0 , while the crude estimates above 
suggested that 10* Mq fllaments might be more likely. It 
should be noted though, that the quantity of cold gas, in 
the simulated tail, is more than capable of producing the 
observed optical luminosity. The simulations also show a lot 
of cold, dense gas in the vicinity of the cloud itself. 

We also highlight the possibility that momentum trans- 
fer between the ICM and cold clouds could be a dynamically 
important process in galaxy clusters. This process could ef- 
fectively couple the energy and momentum of an outflow- 
ing wind to the cold material, thereby dissipating some of 
the energy injected by a central AGN. The presence of cold 
gas may therefore provide an additional feedback mechanism 
that traps energy in the central regions of clusters where it 
is most needed. 
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